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Abstract 

We show explicitly how to construct the quantum Lax pair for systems with 
open boundary conditions. We demonstrate the method by applying it to the 
Heisenberg XXZ model with the general integrable boundary terms. 
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Ever since the introduction of the quantum Lax pair into statistical mechanics by 
McCoy and Wu [jl] and of the notion of a one-parameter family of commuting transfer 
matrices by Baxter 0, a great deal of effort has been expended in the search for spin 
chains with Hamiltonian TC and two-dimensional statistical models with transfer matrix 
r which have the integrability property. 

The quantum inverse scattering method places the theory of completely integrable 
quantum system and solvable statistical models in a unified framework ||. The basis 
to apply it to a completely integrable system is to associate an operator version of an 
auxiliary linear problem Q: 

^ n+1 = L n ^ n , ^ n = AJ> n , (1) 

where L n and A n are matrix operators depending on the spectral parameter u, and a 
dot signifies a time derivative. The consistency condition for equations ([[]) with u— 
yields the Lax pair equation: 

L n = A n+ iL n — L n A n . (2) 

All the solved integrable models appear to imply that a model is completely integrable 
if we can find a Lax pair {L n , A n } such that the Lax equation (0) is equivalent to the 
equation of motion of the model 

Starting with a local Lax operator L n (a matrix acting in the auxiliary space V 
with elements acting in the quantum space h n , at site n), for most quantum integrable 
system the product direct of two Lax operators L n , with different spectral parameters 
satisfy a similarity relation 

1Z(u - v)L n (u) ® L n (v) = L n (v) <g> L n (u)TZ(u - v) (3) 

with JZ a c-number matrix acting in V <8> V. The above relation is referred to as the 
local Yang-Baxter relation and can be represented graphically by figure 1. 




figure 1 
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In terms of the operators L n , the monodromy matrix for a chain with length iV is 
expressed as 

T{u) = L N {u)L N -i{u) ■ ■ ■ L a (u)Li(u) (4) 
and if L' n s with different n commute, we further have 

K{u - v)T{u) ® T(v) = T{v) g> T{u)n{u - v) (5) 

called global Yang-Baxter relation. 

We assume the auxiliary space V and the quantum space h are the same and regard 
the elements of L and 1Z as those of the .R-matrix (TZ = VR with V the permutation 
operator which interchanges spaces of V n with V n -\ ), in terms of which the local 
Yang-Baxter relation (§) has the form 

R 12 {u - v)R 13 {u)R 23 {v) = R 23 {v)R 13 {u)R 12 {u - v), (6) 

where Rlj is a matrix in V x V x V, which acts non-trivially in the spaces Vi and Vj 
only. 

In two-dimensional statistical mechanics, the elements of R may be considered as 
vertices of a vertex model. Baxter first noticed the importance of the relation (^|) in 
that context and regarded it as the solvability condition of the vertex model @: we 
introduce the transfer matrix t(u) as a trace of the monodromy matrix t(u) = trT(w). 
The relation (|^) indicates that there exists a family of commuting transfer matrices 
and that -u-expansion of r(u) gives a set of conserved quantities which are involutive. 

The integrability conditions for a system with open boundary condition are for- 
mulated in order that both the Yang-Baxter equation and the boundary Yang-Baxter 
equations (or reflection equations) are satisfied ||. To a quantum system on a fi- 
nite interval with independent boundary conditions at each end, we have to introduce 
reflection matrices K T (u) to describe such boundary conditions. For a PT-invariant R- 
matrix, the fundamental reflection-factorization relations obeyed by K~(u) and K + (u) 
are 0: 

R 12 {u - v)K{{u)R 21 {u + v)K~{v) = K 2 {v)R 12 {u + v)K^(u)R 21 (u - v) (7) 
which is represented graphically by figure 2 
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figure 2 



For the left boundary we have a similar equation: 

R X2 {-u + v)Kf{u) tl M± x R 2l {-u -v- 2p)M 1 K+(v) t2 
= K+(v) t2 M 1 R 12 (-u-v-2p)Mr 1 K+(u) tl R 21 (-u + v) (8) 

Here U denotes the transposition in the % space; = K± ® 1, Kf = 1 ® K± , etc. 
p is a crossing parameter and M is a crossing matrix, both being specific to a given 
i?-matrix. These relations correspond to the constraint of factorized scattering in the 
presence of a wall. 

There is an automorphism |], |Bj between the K~ and K + namely, given a solution 
K~(u), then K~(—u — pfM is a solution K + {u). 

Given an R -matrix and iiT-matrices satisfying (^) and the corresponding open 
chain transfer matrix t(u) is given by 

t{u) = tr (K + {u)T(u)K-(u)T-\-u)) (9) 

which constitutes a one-parameter commutative family. 

Hamiltonians with boundary terms are obtained from the first derivative of the 
transfer matrix (|9[) |J: 

« V-w 1 AT'(0) tr„[/C(0)g«J . . 

where the local sum is given by 



N d 
V # n ,„_i = a — lnr(«) 



+ const I ill) 

u=0 



The index a denotes the auxiliary space and a is a constant. 
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The corresponding Lax pair formulation may be obtained by direct consideration of 
the equations of motion which can be write in the Lax form considering the following 
operator version of an auxiliary linear problem 



n+l 



J -'n * n 



71 



1,2,..., N 



n = 2,3,..., TV 



(12) 



and 



*i=Qi*i , ^n+i=Qn^n+i- (13) 

Note that Qi(u) and Qn(u) are responsible for the boundary terms in the equations of 
motion of the system. The consistency conditions for these equations are the following 
Lax equations: 

L n = A n+l L n - L n A n for n — 2, ...,7V— 1. (14) 

and 

L\= A 2 Li - L\Q\ and Ln= QnL n - L N A N (15) 

These equations specify the evolution only of L n . It means that A n , Qi and Qn 
have to be determined in terms of L n and H n n _\. 

In this note we will show how to find the Lax pair from the boundary Yang-Baxter 
equations. In this way we are including the general integrable boundary terms in a 
previous work of Zhang where the periodic case was considered. We also apply the 
method to the open XXZ model and compare the result with previous calculations. 

To do this we consider a application of the local Yang-Baxter relation, where the 
vertices are being commuted horizontally rather than vertically: 



7^ n>n _i(e)L n (u + e)L n _i(tt) = L„(m)L„_i(u + e)7£ ni „_i(e) 



(16) 



Here the spectral parameters of the local L-operators differ only by an infinitesimal 
amount e. This relation can be represented graphically by figure 3 




i„-l(»+£) 




*i„,„_i(e) 



figure 3 

This is a form of the star-triangle equation used in 



I. This formulation follows 



also from the non-local quantum Lax pair [|I(], [XT 
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Using the normalization 1Z(0) = 1 we have the following e-expansions 



U{e) ~ 1 + e{oT l H + (31) + o(e 2 
L n {u + e) ~ L n (u) + eL;(w) + o(e 2 ), 



(17) 
(18) 



where a prime stands for w-derivative, with a and (3 some constants and H n n _i of the 
local sum in ([□]) is the .ff in fll7D acting on the quantum spaces at sites n and n — 1. 

Substituting ( |17"D and ([18|) into ([[(]) we get as the first non-trivial consequence the 
following commutation relation 



a 1 [H n , n _ 1 ,L n (u)L n _ 1 (u)] = L n {u)L' n _i(u) - L' n {u)L n _ x {u) 

n = 2,3, ...,N 

We multiply out the above equation to arrive at 

[H n>n -i, L n ] = a(L n L n _ 1 L n ^ 1 — L n ) — L n [H n ^i, L n ]L n \ 
[H n +i, n , L n ] = a(L' n — L n ^_ 1 L' n+1 L n ) — L n+1 [H n +i >n , L n+ i]L n 



(19) 



(20) 
(21) 



forn = 2,3,...,iV-l. 

In order to include the boundary terms into ( p0|) and (|2T|) we also shall consider a 
reflection equation where the vertices are being commuted horizontally: 



L 1 (u)L (u + 2e)K^(e) 



Kz(u + e) 



K^{e)L l (u + 2e)U{u) K^{u + e) (22) 



This equation can be represented graphically by figure 4 




u-\-e 




K~ (e) 



figure 4 

Using the the normalization K~(0) = 1 and 

K~(e) ~ 1 + eK-'(0) + o(e 2 ) 
we get from (p2|) a new commutation relation 



(23) 



-iTr'(O), L a («)L («)] = ^i(«)Lo(«) - L»L ( 



u 



(24) 



We multiply out the above equation and arrive at 

[^-'(OUi] = OWo 1 - L[) - L^K^iO),^ 1 (25) 

Note that in these equations we have extended the chain in order to include the site 
n = at the right boundary and made use of the non-singular property of the L- 
operators. 

Now we recall the equation ( PDj ) to see that the identification 

a~ l H lfi = ^r'(O) (26) 

allows its continuation for n = 1. 

By similar considerations one can derive the left boundary relations. Now the 
equation fl2ip is continued for n = N with the identification 

-Irr tT[K+{G)H Na ] 

a Hn+1 > N ~ trAT+(0) (27) 

These results tell us that we can keep all considerations made by Zhang for the 
periodic caseQ. The main difference consist in remove the term Hn } n+i = H^,i from 
the periodic Hamiltonian and adding two boundary terms H\p and H N+ i N determined 
by the matrices K ± (u). 

The equation of motion for L n is the Heisenberg equation 

L n = i[H, L n ] = i[H n+ i >n , L n ] + i[H ntn -i, Ln}, (28) 

here we have set % — 1. The second equality is because of the locality of the 7i and 
L n . Combining (|2Q|),(|21f) and (p8p, we have 

L n = — ' l L n ] rl \aL' n+1 + [H n +i, n , L n+ ^L n + 'iL n \aL' n _ l — [if n) „_i, L n _i]|L n ^ 1 (29) 



By comparing with the Lax equation flU]), we can read off the second Lax operator 



A n = -LL" 1 {aLj, + [# n , n _i, L n } (30) 



or 



A n = -i {aL'^-lH^, L n ^L-\ (31) 

The compatibility of fl50|) and ( pT| ) is guaranteed by the commutation relation (|20|). 

Therefore, the second Lax operator for completely integrable open chains has the 
following form: In the bulk it is the same for the corresponding periodic chain 

A n = i-f^Ti,n — 1 — iL^Hn n-lLn — iaL n L' n 

n = 2,3,..., N, (32) 
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at the right boundary it is given by 

Qi = ia {\hr(0) - h^K^'iO)^ - L?L' X 
and at the left boundary it is read off from the equation (fill): 

Qn = Mn+i,n — iL N H N+ i N L x 1 — iaL' N L N X 



(33) 



(34) 



where Hn+i,n is given by fl27j). 

Finally, we shall apply this method to a concrete model, the one- dimensional Heisen- 
berg XXZ open chain with Hamiltonian ||, : 



N 



n 



k=2 V Z 

+ sin 2r) (A_al + B_a^ + C_<7^ + + + 



(7 



N 



where 



A, 



cot(^) , E T 



(35) 



(36) 



sin£ T ' sin^ 

Here a x ,a v ,a z and cr^ = ^(cr 1 ± icr y ) are the usual Pauli spin-^ operators, 7/ is a 
parameter associated with the anisotropy, b T , c T and £ T are some constants describing 
the boundary effects. 

It is not difficult to verify that the equations of motion derived from the Hamiltonian 
(j35j) may be cast in the Lax form fll4|) and fll5). 



The first Lax operator L n is identified with the f?-matrix of the 6-vertex model f!4 

/ W4 + w 3 a z n 



R„ 



2w x o n 



\ 2wia+ - w z a z n 

= w 4 + w 3 a z a z n + 2w 1 (a~a+ + a + a-) 

where the elements are parametrized by trigonometric functions of u 

+ 1^3 = sin(u + 27/), u>4 — w-s = sin u, 2w\ = sin 2r\ 

Using the e-expansion ( |T7| ) we have 

a = — sin 2rj and (3 = — - cot 2t/ 
The inverse and the first derivative of L n are well-defined: 



(37) 



(38) 



(39) 



L- 1 = v A + v 3 a z a z n + 2v l (a-a+ + a + a-) 



w' 4 + w' 3 a z a z n 



(40) 
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where 



1 1 

— sm(u-2r)), v 4 -v 3 = —smu, 



v 4 + v 3 = 

A = sin(w — 2rf) sin(-u + 2rf) 



2v\ = — — sin 27? 



(41) 



Now, by direct substitution of these data in ( |3"2] ) one can easily find the second Lax 
operator A n for n = 2, 3, N : 



A r , 



iconst.J + id^Xcr^! + i/(u) (cr n + cr+a^ 



+icr 

-i<7- 

— icr" 



K ff tl- ff nVl 



_p(u) - + q(u) (cr~ + 

_p(u) (<<r+_i - 0-+<-i) + (<r+ + 



(42) 



where 

d(u) 

g(u) 

const. 



1 1 

sin Ar] sin 2r/, f(u) = — sin 77 (sin r\ cos 2u + sin 3u) , 

1 . . 1 . . 1 . t 

— — sin 2ri sin 2u, = — — sin 47? sin u, qiu) = — -sin4?7sinw 

2A ' v ; 2A ' v ; 2A ' 



2A 



sin 2?7(sin 2u — sin 2i] cos 277) 



(43) 



This solution is the Lax operator for the periodic XXZ spin chain which was obtained 
by Sogo and Wadati [15] in the trigonometric limit of the Lax pair operators for the 
one-dimensional XYZ Heisenberg spin chain. 

Now we have to compute the corresponding boundary operators. By a direct cal- 
culation with the operators L x and L^ 1 given by Q4"DD and together with one gets 
the following results for each term of Q\. 



iH 
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— \otL 1 L\ 



-ia (A_crf + B_o\ + C-O v 



— ia {v4w' 4 + — ia (v4W 3 + 1*3 if^) aV^ 

— 2iv\ (w' 4 — w' 3 ) (v~<7i + cr + o"r) 



(44) 



(45) 



and 



—iL^ 1 HqiLi = iaA-^v^Wi + v 3 w 3 — 2v\W\)a\ + iaA^^v^ws + v 3 w± + 2v\W\)a z 
+2iaA_ [(i> 4 — 113)101 — Vi(w4 — w 3 )] (<J~a± — o" + crf ) 
+ia(v^W4 — v 3 w 3 )(B-af + C_af ) + ia(v 3 W4, — V4W 3 )a z (B_af — C_af ) 
+ia [u7i(w 4 + ^3) + vi{w± + w 3 )] (S_cr + + C-o~) 

+ia [wi(v4, + «s) - «i(t«4 + Wa)] (£_ct + - C_a~)<7i (46) 
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We thus find the following Lax operator 
Qi = iconst.l + 



i sin 2 2rj 
A sin 



( — \ sin (2rj + £_)cr^ + \ cos £_ sin 2r\ sin ( u — ^_)cr 1 \ 

+6_t(M)<Ti" + c_r(M)cr 1 " +6_ sin 2rj cos w — c_ sin w cos 2r/crj : 

— sin(n + C-) a i ~ \ sin(2?7 — £_)cr| — | cos£_ sin 2r] 

^ +&_ sin2r]cosii + c_ sin m cos 2r/o"f +6_r(-u)a 1 f + c_t(w)crf y 

(47) 



where 



, sin it cos(w - 77) , sm w cosfw + 77) . _ 

= sin 277 and r(w) = sinir/. (4c 

cos 77 cos 77 



The left boundary operator Q N is obtained using the equation fl3"4"|). This results that 
Qn is the transposition of Qx, followed by the following substitution: 

£_ — > £ + , 6_ — > c + , c_ — > 6 + and o\ — > <7n (49) 

These Lax operators are the trigonometric limit of the Lax pair for the open XYZ spin 
chain given in ref.|T3|. In particular, the diagonal case (6 T = c T = 0), was first derived 
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TBI 



To summarize: solving the Yang-Baxter equation together with the reflection equa- 
tions, we can read off L n , 1Z and K T operators; if Ti has the form (ID) with its bulk 



([[]) is related to 1Z by (|TT|), then the corresponding A n operator is given by (|S2]). More- 
over, for every solution K~ (K + ) of the reflection equation (|7|) (([S])) with K~(0) ^ 
(TrA + (0) ^ 0), we can read off the right (left) boundary term of TC which has the form 
(Eff) ((E7p), then the corresponding Qi (Qn) operator is given by ( |3"B| ) ((0)). 
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